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Abstract. A detailed investigation of the reaction np — » ppn~ has been carried out using the data obtained 
with the continuous neutron beam produced by charge exchange scattering of protons off a deuterium 
target. A partial wave event-by-event based maximum likelihood analysis was applied to determine contri- 
butions of different partial waves to the pion production process. The combined analysis of the np — > ppir~ 
and pp — > ppn° data measured in the same energy region allows us to determine the contribution of isoscalar 
partial waves (1=0) in the momentum range from 1.1 up to 1.8 GeV/c. The decay of isoscalar partial waves 
into CSo) ppTi channel provides a good tool for a determination of the pp S-wave scalar scattering length 
in the final state which was found to be a pp = —7.5 ± 0.3 fm. 

PACS. 13.75.Cs Nucleon-nucleon interactions - 13.85.Lg Total cross sections - 25.40.Fq Inelastic neutron 
scattering 



1 Introduction 

The single pion production in NN collisions is the main 
inelastic process at energies below 1 GeV. Despite the fact 
that a lot of experiments have been performed in this re- 
gion, many questions on this process are still far from be- 
ing satisfactory answered. One of them is the question 
' about contributions of isoscalar (1=0) partial waves to 
, the inelastic neutron-proton collision. The neutron-proton 
scattering amplitude contains both isoscalar and isovector 
(1=0,1) parts, and while the isovector part is rather well 
known, even the magnitude of the total isoscalar cross- 
section is badly determined. Usually, this cross section is 
extracted from the difference of the total cross-sections 
of the pion production reactions: np — ► ppir~ and pp — » 
ppn" . However the numbers obtained in different experi- 
ments are fairly scattered and do not give the whole pic- 
ture for the behavior of the isoscalar cross-section in the 
energy region below 1 GeV. The only conclusion which 
can be made from these values is that the contribution of 
the isoscalar cross section to the np — > ppir~ reaction is 
smaller by about one order of magnitude than the contri- 
bution of the isovector one. 

It should be noted that the experimental data on the 
np — > ppm~ reaction below the energy of 1 Gev are much 
poorer than those for pp collisions. The main reason for 
this situation is due to the difficulty to create a pure mo- 
noenergetic neutron beam. Therefore most previous exper- 
iments used either a continuous neutron beam in a fairly 
wide energy range [1] or investigated the pn-interaction 



with a proton beam and a deuteron target as a neutron 
source [2,3,4 . Of course, it is worth to mention the works 
where the energies of neutrons were defined using the time- 
of- flight technique [5J[B]. 

In the present study we use the continuous neutron 
beam produced by the pd interactions on the deuterium 
target. For a determination of the contributions of the 
isoscalar cross-section we performed a partial wave anal- 
ysis of the data in whole region of incident neutron mo- 
menta combined with the pp — > ppn data. In such an 
approach, the contribution of the isoscalar channel can be 
determined, for example, from the asymmetry of the pion 
distribution in the c.m.s. of the reaction, which is defined 
by the interference of isoscalar and isovector amplitudes. 
Thus the small isoscalar amplitudes can be defined with 
good accuracy. 



2 Experimental set up and data selection 

The np — > ppn~ data were taken at the PNPI synchrocy- 
clotron with the 35 cm hydrogen bubble chamber disposed 
in the 1.48 T magnetic field. The neutrons were produced 
from collision of 1 GeV energy proton beam with thin- 
walled liquid deuterium target. The charged particles are 
swept away from the created neutron beam by a system 
of the additional magnets and lenses situated behind the 
deuteron target. The neutron beam passed through a col- 
limator with size 100x2x1 cm 3 and then the distance of 
5m to the bubble chamber. The produced neutrons had a 
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momentum spread in a range around 1 -j- 2 GeV/c. The 
detailed description of the neutron beam and its energy 
distribution can be found in [7]. 

A total of 10 5 stereofilms were obtained with this beam. 
The films were scanned twice with a goal to select events 
with two positive curvature tracks and one negative track. 
Due to small density of tracks on frames the efficiency of 
such selection was near to 100%. As the result, 10835 three 
prong events were found. Those events can belong either 
to the single pion production reaction 



np — ► ppn 

or to the double pion production reactions 

np — > npir + Tr~ , 
np — > dTr + n~ . 



(1) 



(2) 
(3) 



To separate the first process from the other two, tracks 
of events disposed in the fiducial volume of the chamber 
were measured, geometrically reconstructed and the mo- 
menta of assumed particles were determined in according 
with their masses. The identification of the events was per- 
formed by a kinematic fit demanding that the confidence 
level of accepted events is better then 1%. The direction 
of the incident neutron beam is known from the collimat- 
ing system which provides the direction with the angular 
spread of 0.3° [7]. But nothing is known about a value of 
the neutron momentum which had to be defined by a Semi- 
constrained fit. Three prong events were constrained by 
the kinematical fits for all three mentioned above reactions 
(l)-(3). A visual scan of a bubble density on the tracks 
was carried out for events fitted by several hypotheses in 
order to distinguish between pion, proton and deuteron. 
The angular resolution of the direction of the incident neu- 
tron beam both in azimuthal and in polar angles was first 
estimated to the 1° which turned out to be more than 
three times larger than that obtained after the fit. After 
repeated fits with corrected errors we selected 8251 events 
assigned to reaction (1). In the present experiment we did 
not monitor the neutron beam and therefore the momen- 
tum distribution is given not in milibarns but in the num- 
ber of events. To determine the absolute isoscalar cross 
section we included in our partial wave analysis the pp — > 
ppn data for which the total cross section was measured 
in our previous experiments. 

Fig.[T]shows the neutron beam distribution for selected 
np — > ppn~ events. The neutron momentum distribution 
shows a peak approximately at the energy of the origi- 
nal protons incident on the deuterium target. The second 
peak with less energetic neutrons is created in inelastic 
processes. 

Angular distributions of pions and protons in the c.m.s. 
of the reaction for chosen intervals of the initial neutron 
momenta are shown in Fig.O The proton angular distribu- 
tions are rather symmetrical due to permutation of pro- 
tons, while the pion angular distributions show a clear 
asymmetry (especially at low momenta of incident neu- 
trons). While isovector amplitudes are symmetrical with 
respect to the direction of the produced pion, the isoscalar 
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Fig. 1. The neutron beam distribution for events of the reac- 
tion np — > ppn~ . 

amplitudes are antisymmetrical ones. Therefore, the asym- 
metry in the differential cross section is defined by the in- 
terference of these amplitudes providing a sensitive tool 
for a determination of the weak isoscalar partial waves. 

The two body invariant mass distributions of final par- 
ticles are shown in Fig. [3] The irp invariant mass at high 
energies of the neutrons has a peak defined by the produc- 
tion of Z\(1232) state and indeed, as the partial wave ana- 
lysis has shown, the Z\(1232)p channel is one of dominant 
final states in this reaction. The pp invariant mass distri- 
butions are smooth and do not reveal any peculiarities. 



3 Parameterization of partial waves 

For the production of three particles with the 4-momenta 
qi from two particles colliding with 4-momenta k± and k 2 , 
the cross section is given by: 

do= d<P 3 (P,q u q 2 ,q 3 ) , P = k 1 + k 2 , (4) 

where A is the reaction amplitude, k is the 3-momentum 
of the initial particle calculated in the c.m.s system of the 
reaction, and s — P 2 — (ki + k 2 ) 2 ''. The invariant three 
particle phase space is given by 



d$ m (P, q m ) = 8 i (P - J2 Qi)t[ 



<Pq l 



\W2q 0l 



(5) 



The total amplitude can be written as sum of the partial 
wave amplitudes: 

A = Y J At{s)Q^ J (SLJ)A 2b {i,S 2 L 2 J 2 ){ Sl ) x 

a. 

QiZ.^S 2 L 2 J 2 S'L'J) . (6) 
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Fig. 2. The angular distribution of pion (left-hand column) 
and final proton (right-hand column) in the c.m.s. of the 
np — > ppir~ reaction for 4 intervals of neutron momentum 
(in GeV/c). The full curves are the result of the partial wave 
analysis, the dashed curves are contributions from isovector 
amplitudes and dotted curves are contributions from isoscalar 
amplitudes. 



Here S,L, J are spin, orbital momentum and total angu- 
lar momentum of the np system, S 2 ,L 2l J 2 are spin, or- 
bital momentum and total angular momentum of a two- 
particle system in the final state and S',L' are spin and 
orbital momentum between the two particle system and a 
third particle with momentum qi. The invariant mass of 
the two body system can be calculated as Si = (P — qi) 2 . 
The multiindex a denotes all possible combinations of the 
S, L, J, S 2 , L 2 , J 2 , S', L' and i, A? r (s) is the transition am- 
plitude and A 2 b(i, S 2 L 2 J 2 )(si) describes rescattering pro- 
cesses in the final two particle channel (e.g. production of 
Z\(1232)). In this spin-orbital momentum decomposition 
we follow the formalism given in [8,9,10]. The exact form 
of the operators for initial states QJ 1 "... A1J (S'iJ) and final 
states Q£j% 7 (i, S 2 L 2 J 2 S' ' V 'J) can be found in [TDj. 

Following this decomposition we use the usual spectro- 
scopic notation 2S+1 Lj for description of an initial state, 
the system of two final particles and the system "specta- 
tor and two-final particle state" . For the initial np system, 
states with total momenta J < 2 and angular momenta 
between two nucleons of L — 0,1,2,3 were taken into ac- 
count. For the final three particle system we restricted the 
fit by angular momenta L 2 — 0, 1, 2 and V — 0, 1, 2. 




Fig. 3. The up and pp invariant mass distributions of the 
final particles for 4 intervals of neutron momentum (in GeV/c). 
The full curves are the result of the partial wave analysis, the 
dashed curves are contributions from isovector amplitudes and 
dotted curves contributions from isoscalar amplitudes. 



Due to the nonresonant nature of the np system in the 
energy region investigated here, there is no factorization 
between initial and final vertices and the transition am- 
plitude can depend on all quantum numbers which char- 
acterize a partial wave (index a). Moreover, due to con- 
tribution of the triangle singularities, the production pa- 
rameters can be complex numbers. The best description 
was obtained with the parameterization 



(7) 



where af are real parameters. The a" parameter corre- 
sponds to a pole situated in the region of left-hand side 
singularities of the partial wave amplitudes. It is intro- 
duced to suppress the growing of the amplitudes at the 
regions of large s. 

We also used other, more complicated parameteriza- 
tions of the transition amplitude. However, cither we ob- 
tained a worse description of the data or a similar descrip- 
tion with larger number of parameters in the fit. In the 
latter case we included those results to determine the fi- 
nal systematical errors in the isoscalar contribution to the 
cross section. 

For the irp system in the intermediate state we intro- 
duce two resonances, Z\(1232)i-33 and Roper 7V(1440)Pn. 
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Fig. 4. The angular distribution of the final pion and proton 
in the helicity and Godfrey- Jackson system for the total inter- 
val of neutron momentum. The full histograms are the result of 
the partial wave analysis, the dashed histograms are the contri- 
butions from isovector amplitudes and dotted histograms are 
contributions from isoscalar amplitudes. 



For the 4(1232), we use a relativistic Breit-Wigner for- 
mula with mass and width fixed from the PDG values. 
The Roper state was parameterized in agreement with 
Brcit-Wigner couplings found in the analysis [11]. Let us 
note that the present analysis is not sensitive to the ex- 
act parametrization of the Roper resonance: only the low 
energy tail of this state can influence the data. 

For the description of the final pp interaction we use a 
modified scattering length approximation formula: 



\rPq*4 P + iqa§ p q 2L /F(q, r\ L) ' 



1 



(8) 



where multiindex (3 denotes possible combinations of a 
kinematical channel i and quantum numbers Sa, L2 and 
J 2, dpp is a pp-scattering length and is the effective 
range of the pp system. The F(q,r,L) is the Blatt-Weiss- 
kopf form factor (it is equal to 1 for L = and the explicit 
form for other partial waves can be found, for example in 
[8]) and q is a relative momentum in the pp-system: 



4m^ 



(9) 



For the S-wave partial wave this formula corresponds 
exactly to the scattering length approximation suggested 
in [12 , 13J . To check a possible additional energy depen- 
dence we also introduced a more complicated parameter- 
ization of the numerator in Eq. (jSJ) by substituting: 




1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 
p(n), GeV/c 

Fig. 5. The isoscalar total cross section. The band represents 
the systematical errors obtained from the large number of fits. 
The experimental points are taken from |15] (squares) and [17] 
(triangles) . 

where b 13 are fit parameters. However we found that this 
freedom does not improve the description but leads to 
large correlations between these parameters and parame- 
ters of the transition amplitude Af r (s). 



4 The results and discussion 

We minimized the log-likelihood value fitting simultane- 
ously the present data taken in the whole range of the neu- 
tron momenta and data obtained earlier on pp — ► ppir [14l 
115] measured at nine energies covering the same energy in- 
terval as the np data. To fix the low energy region one set 
of high statistics data taken by the Tubingen group [16] 
on pp — > ppir was also introduced to the fit. The pp colli- 
sion events are produced solely from the isovector channel 
(1 = 1) and therefore the contribution of isovector states 
to the np — ► ppir~ reaction (which differs only by a com- 
mon Clebsch-Gordan coefficient) is strongly constrained. 
This appreciably increases the accuracy for the extraction 
of the isoscalar contribution from the data. 

The results of the partial wave analysis are shown in 
Figs. [2]{3] by solid curves. Contributions of the isovector 
channel are given by dashed curves and contributions from 
isoscalar channel as dotted curves. It is seen that the fit re- 
produces well the asymmetry of the pion angular distribu- 
tion in the overall c.m.s. (see Fig. [2]). Moreover, the angu- 
lar distribution of isoscalar channel changes notably with 
energy demonstrating contributions from different partial 
waves. 

The largest contributions from isovector initial states 
to the pion production process stem from the 3 P2, 3 Pi 
(10) waves (see Table [TJ. At large neutron momenta the 3 J 7 2 
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Table 

( 2S+1 L 



1. Largest contributions of initial partial waves 
j) to the reaction np — > ppir~ . The central values are 



given for the best fit and errors are extracted from the large 
number of different solutions. 
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contribution is also significant. The 1 Sq wave gives the ne- 
gligible contribution. The 3 P 2 state decays almost equal- 
ly into the Z\(1232)p, N(lU0)Pnp and ( 3 P2) PP n states. 
The 3 Pi initial state decays also almost equally into 
Z\(1232)p and ( 3 p2)pp 7r j while 3 i<2 decays dominantly into 
the Z\(1232)p and iV(1440)Pnp final states. The dominant 
isoscalar state 1 Pi decays dominantly into ( 3 Pi) pp 7T and 
jV(1440)Pnp, while 3 Si and 3 D\ states decay almost to- 
tally into the ( 1 So) pp tt channel. 

The quality of the partial wave analysis is also demon- 
strated in Fig. |4l Here, angular distributions of the pion 
and final proton are compared to the fit in the helicity and 
Godfrey- Jackson frames for the total momentum range of 
the incident neutron. The helicity frame is a rest frame of 
the two final particles with the angle calculated between 
one of constituent particles and a spectator particle. This 
frame is most suitable for the investigation of processes 
when two initial particles form a system (e.g. resonance) 
which decays into a final two-body system and a specta- 
tor. The Godfrey-Jackson frame is defined as rest frame 
of two final particles, however here the angle is calculated 
between one of constituent particles and the beam. This 
system is mostly suitable for the study of two particle sys- 
tems in the final state produced by i-channel exchange. 
We would like to mention that the description of these 
distributions at the same momentum intervals as given in 
Figs. I2|3l has the same quality as a total distribution. 

The partial wave analysis defines the relative contri- 
bution of the isovector and isoscalar channels to the total 
cross section. The isovector total cross section was mea- 
sured experimentally [TUfTS] which provides the possibility 
to define the energy dependence of the total cross section 
for the np — > ppix~ reaction. 

A comparison of the isoscalar cross section found in the 
partial wave analysis with results from other experiments 
extracted as: 

a(I = 0) = 3[2<j(np — > ppir~) — <r(pp — > ppn )] , (11) 

is shown in Fig. [5] The points inside the band correspond 
to our best solution. The band itself shows the system- 
atical errors obtained from a large set of solutions with 
different parameterizations of partial wave amplitudes. It 
is seen that our result is systematically exceeds the previ- 
ous calculations. 

The isoscalar cross section increases smoothly from 1.1 
GeV/c up to 1.6 GeV/c where it reaches a value of 2.5 mb 



o 

•*—> 
C 

CD 




cos(7i)h in np 



Fig. 6. The angular distribution of the final pion in the he- 
licity system for two intervals of incident neutron momentum 
(in GeV/c). The full histograms (close to the data points) are 
the result of the partial wave analysis, the dashed histograms 
show the contributions from isovector channel and dotted his- 
tograms the contributions from isoscalar channel decaying into 



the ( 1 S'o)pp7r final state. 



and then it drops sharply down to 1.5 mb. It is possible 
that this drop is connected with the opening of double 
pion production channels where isoscalar np states could 
decay, for example, into two A states. 

The isoscalar initial channel provides a very good tool 
for a determination of the scattering length of the final pp 
system in the pion production reactions. In the initial fits 
we fixed scattering length and effective range of the 1 Sq 



final pp state to a p 



-7.83 fm and r = 2.8 fm [IS]. For 



other partial waves the effective ranges were also fixed to 
r = 2.8 fm while parameters a@ p were fitted freely. 

An appreciable fraction of the isoscalar channel decays 
into pion and the pp system in the final 1 5'o state. The 1 5o 
state produces a sharp peak in the pion backward angle 
in the helicity frame. This angle corresponds to the low- 
est pp invariant mass and due to kinematical conditions 
(suppression of other partial waves), this peak is well ob- 
served in the data. The pion angular distributions in the 
helicity frame for two intervals of the incident neutron 
momentum is shown in Fig. [6] It is seen that the isoscalar 
channel produces a dominant effect at small and large mo- 
menta of incident neutron. If the effective range is fixed 
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at 2.8 fm, we found the scattering length to be equal to 
a vv = —7.5 ± 0.3 fm. 



5 Conclusions 

We performed a detailed study of the np — ► ppir~ reac- 
tion. The data were analyzed in a maximum likelihood 
partial wave analysis. The partial wave analysis provides 
a solid determination of weak isoscalar amplitudes from 
their interference with strong isovector amplitudes. The 
dominant decay of the initial isoscalar states 3 Si and 3 Di 
into the ( 1 So) pp tt final state allows us to obtain a good de- 
termination of the pp scattering length in the final state 
of pion production reactions. 
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